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SECTION 


INTRODUCTION 

In  many  practical  problems  in  hypersonic,  flight,  bodies  having 
surface  discontinuities  or  regions  with  rapid  variation  of  curvature 
are  used.  Typical  bodies  of  such  type  are,  for  example,  cone-cylinder 
combinations . 

A  sublayer  model  for  the  boundary  layer  characteristics  down¬ 
stream  of  the.  corner  was  first  introduced  by  Sternberg*.  In  • 

Sternberg's  paper,  the  predicted  boundary  layer  characteristics  are 
compared  with  measurements  of  the  recovery  factor.  .At  a  later  date, 
independently,  Zakkay  and  Tani^  reintroduced  the  same. model  and  pre¬ 
sented  a  detailed  method  for  evaluating  the  development  of  shear  Layer 
and  sublayer  as  well  as  the  heat  transfer  rate  when  the  oncoming 
boundary  layer  is  laminar.  These  analyses  were  made  by  using  an 
expansion  technique  similar  to  that  used  by  GSrtler.  In  the  present 
work,  the  above  technique  is  extended  to  include  both  transitional  and 
turbulent  boundary  layers.  Again.  In  this  case,  the  method  chosen  for 
analyzing  the  boundary  layer  characteristics  is  based  on  the  model 
introduced  In  reference  1. 

According  to  the  sublayer  model,  the  flow  field  of  interest  may 
be  divided  into  two  layers,  an  inner  nonsimilar  layer  starting  at  the 
corner,  and  an  outer  shear  layer.  The  initial  conditions  for  the  shear 
layer  are  obtained  by  the  invlsoid  and  adiabatic  expansion  of  the  boundary 
layer  profiles  ahead  of  the  corner.  The.  new  viscous  layer  is  shown  in 
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Fig,  1  and  will  be  referred  to  as  sublayer  since  it  underlies  the  outer 
shear  layer.  Distinction  should  be  made  between  the  proposed  sublayer 
and  the  one  usually  associated  with  turbulent  flow. 

This  flow  model  presupposes  that  the  initial  velocity  profile  of 

the  shear  layer  is  very  steep.  In  effect,  the  flow  field  downstream  of 

the  corner  will  be  considered  the  same  as  that  which  would  exist  over 

a  flat  plate,  with  oncoming  viscous  flow  having  shear  and  shear  gradient. 

Consequently,  the  corner  itself  is  treated  as  a  singularity  similar  to  the 

leading  edge  of  a  flat  plate.  Hence,  the  problem  is  similar  to  the  one 
3  4 

treated  by  Li  and  Glauert  .  However,  in  this  case,  the  analysis  uses 
the  experimental  pressure  distribution  and  will  not  be  concerned  with 
the  higher  order  effects  such  as  the  induced  pressure  gradient. 

When  the  boundary  layer  flow  is  fully  turbulent  up  to  the  corner, 
the  sublayer  model  is  expected  to  become  more  effective,  since  the 
velocity  profile  is  steeper  than  the  one  obtained  for  laminar  flow. 

Mathematical  treatment,  however,  is  not  developed  rigorously 
because  of  the  new  unknowns,  i.  e.  ,  Reynolds  stress  and  turbulent  heat 
flux.  Moreover,  according  to  the  equivalence  of  the  flow  field  to  a  flow 
over  a  flat  plate  placed  in  a  turbulent  oncoming  flow,  the  flow  repre¬ 
sented  by  the  sublayer  may  be  laminar  close  to  the  corner  and  may 
undergo  the  transition  to  turbulence  somewhere  downstream.  Four 
independent  sets  of  experimental  data  presented  show  that  the  observed 
heat  transfer  rate  has  such  a  trend.  Taking  into  account  such  a 
phenomenon,  the  sublayer  plus  simple  mass  flow  consideration  gives 
good  results  of  heat  transfer  rate  near  the  corner  and  somewhat  down¬ 
stream  of  the  corner  where  the.  state  of  affairs  becomes  asymptotic. 
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By  these  relations  Eq.  (5)  immediately  follows  from  Eq,  (4), 


For  the.  sublayer,  the  stream  function  and  the  enthalpy  are 
expanded  as  follows: 

iH s » n)  =  (2a)"5  [f  (n)  +  pj.fi(n)  +  a  {eiafn(ri)  +  BiMn)  }+•••] 
H(s,  q )  =  hw  +  He[Y0g0(q)  +  b~  {  SiY0fi  (tl)  +  YlgiW  } 

+  8  {sis Y0^ii  (“H)  +  eigi(ri)  }+ •••] 


(7) 


(8) 


where  @i  ,  Bi  ,  .  .  .  ,  V  ,  Yi  ,  Bi  ,  .  .  ,  are  constants  and  H  is  a  refer - 
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ence  enthalpy.  These  expressions  arc  substituted  into  the  governing 
differential  equations  and  the  coefficients  of  the  same  power  of  b  are 
equated  to  zero.  Thus,  one  arrives  at  a  set  of  differential,  equations 
for  i'o,  fi ,  etc.  The  boundary  conditions  at  the  wall  are  the  no-slip 
and  the  constant  temperature  (H  =  hw).  The  outer  boundary  conditions 
are  found  tay  matching  the  profiles  with  those  given  by  Eqs.  (4)  and  (5) 
at  the  edge  of  the  sublayer.  Thus  the  following  sets  of  differential 
equations  are  obtained: 
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is  a  constant  to  be  determined  by  the  behavior  of  fi  such  that 
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The  experiments  for  a  sharp  cone-cylinder  as  well  as  for  a 
blunt  cone  are  described  in  reference  2.  The  pertinent  flow  quantities 
for  the  20  degree  sharp  cone-cyliner  are: 


Stagnation  pressure  .  =  600  psia 

Stagnation  temperature  =  1800°R 

WalL-to-atagnation  temperature  ratio  =  0.313 

Reg  =  3  x  104 

Free  stream  Mach  number  =  8 

Mach  number  outside  the  boundary  =  4,25  before  the  corner 
layer  6,63  after  the  corner 


Using  the  experimental  data  q  is  obtained  from  Eq.  (8).  Then  Eq.  (19) 
finalLy  results  in 


-^r  .  103  =  +  6.6  +  1 .14  y7  +  •  •  •  (20) 

Reg2  */x” 

where  pt  is  assumed  zero  and  x  is  measured  in  inches.  The  resuLts 
are  plotted  in  Fig.  2.  Also  included  for  comparison  is  the  flat  plate 
solution  corresponding  to  the  uniform  flow  at  the  edge  of  the  outer  shear 
layer.  The  ratio  of  the  fourth  term  in  Eq,  (20),  which  is  omitted,  to  the 
third  term  is  O(»/x/10)  and  is  also  proportional  to  Reg”^,  For  the 
present  value  of  Re^,  the  fourth  term  becomes  of  appreciable  magni¬ 
tude  at  x=10Q  inches. 
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Eqs,  (21)  become 
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where  use  has  been  made  of  unit.y  of  the  turbulent  Prandtl  number 


Fr, 


_ _  Sh 

p  u  'v  '  8y 

^  "p  h'  v  '  8u 
8y 


Upon  comparing  these  equations,  the  Crocco  relation  is  obtained. 


H  -h 
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From  which  the  Reynolds  analogy  results  directly. 
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where  c  /2-p;,,  ~  /p*>i  *  • 
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If  the  empirical  correction  is  made  to  Prandtl  number  being 
other  than  unity,  then 
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Footnote  for  page  16. 

After  completion  of  the  manuscript,  the  authors'  attention  has 
been  called  to  a  skin  friction  formula  given  by  Spence1  . 


where 
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It  can  be  shown  that  this  formula  is  exactly  identical  with  Eq.  (27). 

In  the  reference  state,  the  thickness  of  the  boundary  layer  is  given  by 
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Eq.  (27)  is  arrived  at  if  the  above  expression  is  inserted  into  the 
foregoing  formula. 


(End  of  footnote.) 
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(28) 


is  one -half  of  the  skin-friction  coefficient  of  an  incompressible  flow 
with  the  flow  quantities  corresponding  to  the  outer  edge,  conditions,  the 
factor 
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n+T 


is  the  one  connecting  the  compressible  skin  friction  coefficient  to 
the  incompressible  reference  state.  For  n=7,  (tl*/ue)^’^(p^/pe)^'^ 
is  the  factor  usually  used  in  the  reference  enthalpy  method^. 
Accordingly,  the  undetermined  reference  state  is  taken  as  given  by 


h*  =  0.5(h  +h  )  +  0.22(h  -h  )  (29) 

♦  we  '  aw  e  ' 

1  uc» 

with  haw  =  he+  Fra^  — .  Further  support  for  this  step  may  be  found 

in  the  work  of  Burggraf^  who  proceeded  to  substantiate  the  method.  He 

proposes  the  reference  state  as  that  of  the  outer  edge  of  the  laminar 

sublayer  which  was  first  suggested  by  R.ott^. 

For  a  not  so-slender  cone  in  a  supersonic  flow  with  attached 

14 

bow  shock  wave,  Van  Driest's  law  may  be  applied  to  calculate  the  skin 
friction  and  hence  the  heat  transfer.  According  to  this  rule,  the  cone 
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skin  friction  is  evaluat  ed  with  one-half  the.  Reynolds  number  for  fiat 
plate  with  the  same  free  stream  Mach  number  and  the  same  wall -to- 
free-stream  temperature  ratio.  Thus,  from  Eqs.  (27)  and  (28)  one 
gets 


2  cf 


n+3  u^n+3  pA 
=  2  (— )  {— ) 


cone 


4  w. 

2  fmcomp, 

flat  plate 


(30) 


The  correction  factor  in  this  expression  can  be  written  in  terms  of 
the  Mach  number  and  the  wall-to-free  -stream  temperature  ratio. 
For  example,  assuming  rt  =  7  and 


one  obtains 


U* 


T* 

Gy-) 

e 


0.75 


0.2  M*  0.2  P*  0.8 

2  W  W 


B .  Expansion  Around,  the  Corner 

In  this  step  the  viscosity  is  neglected.  The  velocity  and  stag¬ 
nation  enthalpy  profiles  are  expanded  isentropically  around  the  corner. 
3h  this  manner  a  velocity  profile  having  a  non -zero  value  at  the  wall 
and  a  corresponding  stretched -out  stagnation  enthalpy  profile  is 
obtained  downstream  of  the  corner.  This  scheme  has  been  proposed 
in  reference  1,  and  has  been  used  successfully  in  reference  2. 
Reference  .1  points  out  that  the  high  acceleration  present  at  the  corner 
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shear  Layers  after  the  expansion. 

The  concept  given  above  is  a  result  of  dealing  with  the  micro¬ 
scale  turbulence  Level  present  in  the  boundary  layer  before  and  after 
the  corner.  These  results  have  not  been  confirmed  substantially  and 
are  not  amenable  to  be  used  in  calculating  the  heat  transfer.  There¬ 
fore,  in  this  paper  in  a  manner  analogous  to  that  of  determining  tran¬ 
sition  over  a  flate  plate,  the  momentum  thickness  wiLL  be  used  as  a 
measure  of  the  transition  to  turbulence. 

Therefore,  it  is  important  for  this  case  to  calculate  the  Reynolds 
number  based  on  the  momentum  thickness  before  and  after  the  corner. 

It  will  also  be  demonstrated  in  a  subsequent  section  that  Re^  is  indicative 
of  the  extent  of  the  Laminar  sublayer  after  the  corner. 

Therefore,  since  the  initial  part  of  the  boundary  sublayer  after 
the  corner  is  lamina  r.  araivs.s  dc'  (  'oped  tor  the  laminar  sublayer  is 
also  applicable  here.  The  shear  p^ohies  tor  the  ’aminar  sublayer  will 
be  the  ones  obtained  f*cm  ’he-  mvTs>.  id  expansion  ot  the  turbulent:  profiles 
obtained  ahead  ot  the  corner.  In  this  ,.t  is  assumed  that  the  fluctuations 
in  Ihe  shear  layer  ha'  e  a  negligible  effect  on  the  formulation  of  the 
laminar  Sublayer. 

Next,  consideration  will  be  giver;  ‘c  the  development  of  the 
asvmp'otic  flow  downstream  cf  the  corner.  It.  is  well-known  that  turbu¬ 
lent  boundary  layer  profiles  at  h-.gh  Reynolds  number  preserve  their 
self -sirm iar  character  over  the  greater  part  of  the  flow  provided  there 
i?  no  pressure  gradient^'''7.  For  the  present  case,  despite  the 
negligibly  small  pressure  gradient,  self- similarity  may  not:  be  expected 
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initial  profiles  of  the  outer  shear  layer,  the  local  outer  edge  conditions 
of  the  sublayer  are  determined.  From  Eq.  (33)  the  local  skin  friction 
and  hence  the  local  heat  transfer  rate  are  obtained. 


SECTION  IV 

DISCUSSION  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


In  order  to  demonstrate  the  various  theoretical  aspects  discussed 
in  the  previous  section,  four  experimental  results  are  presented  below. 


Case  1. 

Stagnation  pressure 

Stagnation  temperature 

Wall-to-stagnation  temperature 
ratio 


Re  = 
8 


Pse^eT  *v 


(R  =0.582") 


Free  stream  Mach  number 


Mach  number  outside  the 
boundary  layer 


=  600  psia 
=  1700°R 

=  0.294 

=  4.94x  I04 

=  8 

=  1,70  before  the  corner 
2.67  after  the  corner 


The.  above  test  was  conducted  at  the  hypersonic  facility  of  the 
Polytechnic  Institute  of  Brooklyn,  The  body  consists  of  a  spherically 
capped  24°  half  angle  cone-cylinder,  with  a  bluntness  ratio  of 
RQ/r  =  0.317  (Fig.  4).  The  laminar  boundary  layer  study  over  this 
body  has  been  presented  in  reference  2.  Experimental  laminar  heat 
transfer  results,  as  well  as  comparison  with  theory  before  and  after 
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Case  3 


Stagnation  pressure 

Stagnation  temperature 

Wall-to-stagnation  temperature 
ratio 

Re 

a 

Free  stream  Mach  number 


139  psia 
710°R 

0.47 

3.73  x 10B 
5.0 


Mach  number  outside  the 

boundary  Layer  =  3.8  before  the  corner 

5.12  after  the  corner 


The  results  of  these  tests  have  been  taken  from  reference  20. 
The  body  consisted  of  a  15°  half  angle  sharp  cone  followed  by  a 
cylindrical  afterbody  1.8"  in  diameter.  The  Laminar  results  for  this 
configuration  have  been  presented  in  the  Laminar  portion  of  this  paper. 
In  order  to  obtain  a  turbulent  boundary  layer,  trips  were  placed  on  the 
conical  portion  of  the  body.  The  results  of  the  heat  transfer  meas¬ 
urements  are  presented  in  Fig,  7  in  terms  of  the  symbols  used  in  this 
report.  The  Reynolds  number  based  on  the  turbulent  momentum 
thickness  ahead  of  the  corner  was  calculated  to  be  1370,  while  the  one 
calculated  after  the  expansion  was  1530.  The  analysis  for  the  laminar 
calculations  as  well  as  for  the  turbulent  asymptotic  solution  are 

included  in  Fig.  7.  For  this  case  the  Re0  after  the  corner  is  suffi- 

0 

ciently  large  so  that  the  Laminar  sublayer  is  swallowed  in  a  very  short 
distance  after  the  corner.  Therefore,  the  boundary  layer  is  closely 
approximated  by  the  asymptotic  turbulent  solution.  These  consider¬ 
ations  are  clearly  seen  from  the  theoretical  predictions  included  in 
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Fig.  7.  The  experimental  results  are  much  higher  than  the  laminar 
predictions,  and  are  closer  to  the  turbulent  asymptotic  solution. 


Case  4. 

In  order  to  further  substantiate  the  fact  that  the  extend  of  the 
laminar  sublayer  after  the  corner  may  be  obtained  from  Re^, 
the  experimental  results  of  reference  l  are  used. 


Stagnation  pressure 

Stagnation  temperature 

W a  LI  -  to  stagn  at  i  on  t  emp  e  r  at  u  r  e 
ratio 


=••  27.1  psia 
=  535°R 

=  adiabatic.  walL 


Re 

s 


=  5 . 3  5  2  x  1 0e 


Free  stream  Mach  number  =3.02 


Mach  number  outside  the 

boundary  Layer  ~  1.89  before  the  corner 

3.13  after  the  corner 

The  body  considered  in  the  above  test  consisted  of  a  29°  half 

angle  sharp  cone.  The  adiabatic  wall  temperature  was  measured  before 

and  after  the  corner  for  both  laminar  as  well  as  turbulent  conditions 

ahead  of  the  corner.  For  the  case  where  a  laminar  region  prevailed 

ahead  of!  the  corner.  Re.  is  calculated  to  be  540.  Inviscid  profiles  after 

0 

the  corner  were  also  obtained,  and  the  corresponding  Re^  calculated  for 
these  profiles  was  4h0.  It  is  indicated  in  reference  1  that  transition,  for 
this  case  was  observed  to  occur  at  a  distance  of  20  centimeters  down¬ 
stream  of  the  shoulder.  Since  for  this  condition  the  heat  transfer  is 
zero,  the  momentum  thickness  is  calculated  from  the  ekin  friction 
coefficient,  which  is  ir.  turn  calculated  from  the  theoretical  predictions 
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given  in  the  Laminar  section  of  the  paper,  A  value  of  Re^  equal  to 
1200  was  obtained  from  the  analysis  at  the  observed  station  of  tran¬ 
sition. 

For  the  case  where  the  boundary  layer  was  turbulent  ahead  of 
the  corner,  the  value  of  Re^  calculated  before  and  after  the  corner  is 
1200  and  800  respectively.  The  portion  where  transition  was  observed 
for  this  condition  was  at  3  centimeters  downstream  of  the  shoulder, 
the  Re^  calculated  at  the  point  of  transition  was  equal  to  1100.  There¬ 
fore,  for  this  condition  the  boundary  layer  is  predominantly  turbulent 
as  indicated  in  the  Turbulent  Flow  section  of  this  paper.  What  is  more 
important  is  the  fact  that  transition  occurred  approximately  at  the 
same  value  of  Reg  as  for  the  laminar  case.  Therefore,  the  assumption 
of  neglecting  the  fluctuations  in  the  shear  layer  in  analyzing  the 
boundary  layer  and  the  criterion  that  Reg  is  indicative  of  transition 
in  this  case,  has  clearly  been  demonstrated. 

SECTION  V 
CONCLUSIONS 

A  flow  model  has  been  adopted  for  treating  the  boundary  layer 
downstream  of  a  sharp  corner.  The  analysis  based  on  this  model  for 
a  laminar  flow  was  presented  in  reference  2,  Improvements  are  made 
herein  whereby  a  unified  treatment  of  velocity  and  thermal  profiles 
based  on  the  Crocco  relation  hoLds,  thus  are  amenable  to  the  calculation 
of  the  higher  order  solutions  including  the  pressure  gradient  as 
perturbation. 
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An  extension  of  this  flow  model  to  conditions  where  the  boundary 
layer  ahead  of  the  corner  is  turbulent,  is  also  included.  These  results 
were  compared  with  some  of  the  experimental  results  available  in  the 
literature. 

The  results  also  reaffirmed  the  considerations  presented  by 
Sternberg  that  a  new  laminar  boundary  layer  does  start  at  the  discon¬ 
tinuity.  This  boundary  layer  will  prevail  to '  a  distance,  and  a  second 
transition  to  turbulent  will  subsequently  occur. 

It  is  indicated  here  that  for  the  range  of  test  conditions  observed, 
that  if  Re^  before  the  discontinuity  is  of  the  order  of  600  -*700,  then  a 
turbulent  boundary  layer  exists  ahead  of  the  discontinuity.  The  expan¬ 
sion  around  the  corner  has  a  *  wo -fold  die,  j  on  Re  .  The  first  ir  an 

0 

increase  ir.  the  momentum  kne*«.  0,  ar.d  a  decrease  in  the  density. 

Based  on  the  caJci.jat.OTi  of  Re  at*-r  *  ho  v.crner,  ar  idea  where  frar.- 

c 

sit-ion  occurs  on  the  ..yhr.de  r  mav  be  obia  ued.  If  is  irdicated  from  the 

results  presented  here  that  fi.e  se.cr.d  t :  ansitior  ‘o  turbulent  occurred 

at  an  Re  of  the  order  of.  1  L00-.  I  i00.  For  the  case  where  Re  is 

0  9 

larger  than  1400  after  the  expans  or,  the  heat  transfer  is  predominantly 
turbulent  and  may  be  estimated  from  the  asymptotic  turbulent  analysis. 
Based  on  these  ■  or.®  Me -a*  ions,  calculations  for  the  full  laminar,  tran¬ 
sitional-  and  tnrbuierd  boundary  'aver  after  the  discontinuity  are  then 
presented.  !*■  ,?  also  indicated  here  than  the  effect  of  shear  on  the 
development  of  the  turbulent  boundary  taver  after  the  discontinuity  is 
of  much  Jess  significance  than  for  the  laminar  case. 
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TABLE  2 


T) 

gi 

gi 

z 

0 

0 

0.4572 

0.2 

0.0915 

0.4575 

0.4 

0.1831 

0.4595 

0.6 

0.2755 

0,4649 

0.8 

0.3694 

0.4752 

1.0 

0.4660 

0.4918 

1.2 

0.5666 

0.5154 

1.4 

0.6726 

0.5465 

1.6 

0.7856 

0.5846 

1.8 

0.9069 

0.6286 

2.0 

1.0374 

0.6767 

2.2 

1.1777 

0.7267 

2.4 

1.3280 

0.7760 

2.6 

1.4879 

0.8223 

2.8 

1 ,6566 

0.8639 

3.0 

1.8331 

0.8994 

3.2 

2.0159 

0.9284 

3.4 

2.2040 

0.9509 

3.6 

2.3959 

0.9676 

3.8 

2.5907 

0.9795 

4.0 

2.7874 

0.9875 

4.2 

2.9855 

0.9927 

4.4 

3.1844 

0.9959 

4.6 

3.3838 

0.9978 

4.8 

3.5834 

0.9988 

5.0 

3.7833 

0.9994 

5.2 

3.9832 

0.9997 

5.4 

4.1831 

0.9998 

5.6 

4.3831 

0.9999 

5.8 

4.5831 

0.9999 

6.0 

4.7831 

1,0000 
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TABLE  4  (fii  for  p!  =  0) 

i  a 


r\ 

fll 

2  2 

fll' 

2  2 

fn" 

2  2 

0 

0 

0 

0.1576 

0.2 

0.0031 

0.0316 

0.1576 

0,4 

0.0126 

0.0630 

0.1576 

0.6 

0.0283 

0.0946 

0.1578 

0.8 

0.0505 

0,1262 

0.1591 

1.0 

0.0788 

0.1584 

0.1631 

1.2 

0,1138 

0.1919 

0.1730 

1.4 

0.1559 

0.2282 

0.1933 

1.6 

0.2055 

0.2701 

0.2287 

1.8 

0.2645 

0.3210 

0.2829 

2.0 

0.3348 

0.3846 

0.3560 

2.2 

0.4194 

0.4644 

0.4434 

2.4 

0.5217 

0.5623 

0.5350 

2.6 

0.6455 

0.6778 

0.6179 

2.8 

0.7939 

0.8078 

0.6780 

3.0 

0.9692 

0.9468 

0.7047 

3.2 

1.1727 

1.0872 

0.6932 

3.4 

1.4036 

1.2216 

0.6452 

3.6 

1.6604 

1.3433 

0.5687 

3.8 

1.9398 

1.4479 

0.4752 

4.0 

2.2383 

1.5331 

0.3766 

4.2 

2.5518 

1.5989 

0.2836 

4.4 

2.8767 

1.6474 

0.2030 

4.6 

3.2098 

1.6812 

0,1383 

4.8 

3.5485 

1.7038 

0.0898 

5.0 

3.8908 

1.7181 

0.0557 

5.2 

4.2353 

1.7267 

0,0331 

5.4 

4,5813 

1.7319 

0.0191 

5.6 

4.9279 

1.7348 

0,0107 

5.8 

5.2750 

1,7365 

0,0061 

6.0 

5.6224 

1.7374 

0.0036 
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FIG.  4  HEAT  TRANSFER  DISTRIBUTION  FOR  A  24*  BLUNTED  CONE-CYLINDER 


FI 6. 4  a  TURBULENT  BOUNDARY  LAYER  PROFILES 
BEFORE  EXPANSION  FOR  A  24°  BLUNTED 
CONE -CYLINDER. 
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FIG.  7  HEAT  TRANSFER  DISTRIBUTION  FOR 
15°  SHARP  CONE- CYLINDER. 
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U  r»  y 

n  =  — S-t  pdy 

(2e)2  o 


Development  of  the  shear  layer  in  the  neighborhood  of  s  =  0  is 

x 

determined  by  expanding  the  stream  function  in  power  series  of  s2  and 
integrating  the  momentum  and  energy  equations.  The  results, 
expressed  'in  terms  of  ri,  are 


u  =  uo[lh(2s)2  W(r±-  n+  -jf)  +  s  {2a^ 


TT  ki3  . 

+  ki  -5^2-  q  +  —  +  - 

2  ue6i  2  uouP*i' 


i  PiA 

- — r)  +  krlO  - - 

,  u  8i*  p  ,cu  u 

o  e  1  Hei  e  o 


a 

}  +  0(s8)  ] 


H  =  h 


.  Jl  +  (2s)*n<—  n+  -Jr)  +  S  {zi(~~  rf 
0  ue6i  -JZ  k  A,su  3 


k.  3 

.  .  JT  .  i  . 

+  - -  q+  -  + 


2  ue6i 


—  +  - - — r)  +  kiO+  0(sa )  ] 

2  u  u  6i 
o  e 


1  U  3 

where  A  =  —4 —  (h  -  -2-)  and  ki  ,  ki  are  constants.  It  is  observed  that  if 
u  h  o  2  2 

o  ei 

Pi  =  0  one  can  expect  Crocco  relation  between  the  two  profiles  to  hold. 

In  fact,  then  one  can  also  expect  the  same  relation  for  the  initial  pro¬ 
files,  namely,  with  a  constant  Ci  ,  one  gets 

\ 

C,u  +  h  =  h 
1  o  w  o 

CiUouj  =  hQQ  \  •  (6) 

Ciuo°  =  h0i 
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In  terms  of  these  functions,  one  arrives  at  the  following 
formulae,  in  which  some  anomalous  expressions  in  reference  2  have 
been  corrected. 


u  =  u  [f  '  (r)  +  -^5^~  ajf  1  'M  +  — — —  {  ujaf  1 1  /(T',)  +  ofi '  M  }  f  0(sa)  ]  (16) 

e  °  ue4i  2  '  u! 36ls  L  25  J 


H  =  hw+  Ws» w  +  7X  f(ho*hw)«£i/(t,) 


Q 

+  {hon-  (ho-hw) »  }g  1 M  ]+  —~7  {(h0-hw)uj8f|i;  (T1)  +Hee> (r) }  +  °(8a) 


2  V«i‘ 


Note  that  H  is  still  undetermined.  When  pi  is  sufficiently  small, 
e 

one  may  expect  the  Crocco  relation  for  the  initial  profiles,  i,  e,  ,  Eq.  (6). 
This  is  the  case  for  the  profiles  obtained  by  the  inviscid  and  adiabatic 
expansion  technique.  Then  h.Q f3 - (ho -h^) uj  =  0.  Thus  the  last  group  of 
terms  on  the  right-hand  side  of  the  differential  equation  in  Eq.  (15)  dis¬ 


appears.  Choose  such  that 
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H  =  o(h  -h  )  (18) 

e  o  w 

Then  the  differential  equations,  Eq.  (I6')  for  and  Eq.  (13)  for  ^ ' 
become  identical  with  identical  boundary  conditions  provided  pi  =0, 
because  from  Eqs.  (6)  and  {18)  one  has 

H  a( h  -h  ) 

e  _  o  w  _  , 

h  $  "  h  $ 

o  o 

Also,  then  it  can  be  seen  from  Eqs.  (16)  and  (17)  that  the  Crocco 
relation  exists  for  t.he  inner  profiles. 

When  pi  is  not  small,  the  same  relation  cannot  be  expected  for 
the  outer  and  inner  profiles!,  Then  H ^  may  assume  a  value  other  than 
given  by  Eq.  (18). 

Tables  of  the  functions  fi,  gi.,  .  .  .  are  given  at  the  end  of  the 

2  2 

* 

report. 

The  heat,  transfer  rate  at  the  wall  can  be  expressed  in  a  non- 
dimensional  form  as 


Nu 

- X 

Re  2 
s 


qC  R 

P  o  , 

*se _ _ 1 

x  (H  .-h  )  Re  2 
se  ej  w  s 


(19) 


where 


H  ,BH  0n  H 

q  "  c;  W  -ByV-0  - 


uopwUw  ,8H, 

i-  /  >  \y  5"n  n=0 

p^w  (^s)'i 


where  q  is  the  amount  of  heat  transferred  per  unit  lime  area, 
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functions. 
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where  the  boundary  layer  is  artificially  tripped,  and  the  second  case 
will  be  for  the  condition  where  the  free  stream  Reynolds  number  is 
sufficiently  high  for  a  turbulent  boundary  layer  to  be  established  before 
the  corner. 

A.  Method  of  Analyzing  Boundary  Layer  Upstream  of  the  Corner 

In  order  to  insure  that  the  boundary  layer  ahead  of  the  corner 
was  turbulent,  the  heat  transfer  rate  in  each  case  was  calculated  using 
standard  turbulent  boundary  layer  analysis.  The  velocity  and  enthalpy 
profiles  ahead  of  the  corner  were  calculated  by  using  the  tabulations  of 
reference  7  for  a  (l/7)th  power  Law  for  the  velocity  profile. 

In  what  follows,  a  formula  for  the  heat  transfer  rate  on  a  cone, 
Eqs.  (2b)  and  (10)  and  that  for  a  flat  plate,  Eqs.  (26)  through  (28)  are 
derived.  These  formulae  will  be  used  for  the  conical  part  of  the  body 
and  for  the  cylindrical,  afterbody  respectively.  The  main  procedure  is 
such  that  the  flow  field  is  t  ransformed  by  the  Howarth  transformation 
into  a  constant  density  reference  state  which  may  be  regarded  as 
incompressible  flow  field  so  that  the  empirical  power  law  profile  for  a 
flat  plate  may  be  applied.  Also,  the  reference  enthalpy  method  is 
naturally  arrived  at.  For  a  cone.  Van  Driest' a  law  is  applicable  to 
convert  the  flat,  plate,  skin,  friction  values. 

The  boundary  layer  equations  for  a  turbulent  flow  over  a  flat 


plate  are  given  by 


h (pu)  +  ^ (pv)  =  0 


p  U  +  pv  -J—  =  •—  [  (i  -  p  u'v'] 
M  8x  8y  8yl^8y  ^ 


-  9H  .  —  8H  8  r—  8H  -  71 

PuS£-+Pv-^=-gy[|i5^--p  Hv7] 


B  D 

where  ba£  implies  time  average,  ALso  it  is  assumed  that  ^  =  0  and 


cp|J  _ 


PrT  =  -  1.  Introduce  the  stream  function  such  that 

K 


P  U  =  P 


8  1J1 


#  0y 


pv 


Sty 


Sx 


where  star  denotes  the  unspecified  reference  conditions.  Then  the 
Howarth  transformation 

y  7 


y*  s  f  dy 

*  i  p* 


8 


is  applied,  Following  Mager  ,  the  following  quantities  are  assumed 
invariant . 

(1)  the  stream  function  ty  =  ty,.,  where 


St* 


8*,; 


8y* 


v>:<  =  ‘  0X„, 


(2)  the  turbulent  shear  per  elemental  mass 
P  (  u  'V  ')  =  p,((  u  'v  % 


M  =  P 


.  S  Cf 

3  P)|< 

St  =  Pr - r- 

av  PA  2 


(24') 


The  Blasius  skin  friction  formula  for  (1/n)^  power  profile  gives 


-Hrf 


2 

>  nt  3 


2 


{c(n)  ]■ 


2n  ^  (n+  2)  (n+  3) 
n+3 


Re 


(25) 


U  x 

where  Re*  =  — —  and  c(n)  is  given  in  reference  9.  Combining  Eqs.  (25) 
x  ^ 

and  (24),  the  following  expression  for  the  Nueselt  number  is  obtained. 


Nu  =  St  •  Rex  •  Prav  =  Pr&v  .  —  •  Re} 


(26) 


where 


1 

2  Cf=  p  u  » 

b 


2  2  2  1-  2 

iJ+T  nfT  n+5  n+  $ 

1 _ { V-r)  I—)  ft  ft 

2n  l(n+2)(n+3)  J  p  u  x  n@  pe 


(27) 


{c(n)  ]■ 


?T+T 


See  footnote  on  page 
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within  at  least  some  distance  downstream  of  the  corner.  This  is  due 
to  the  transient  state  immediately  after  the  rapid  expansion.  The  flow 
must  adjust  itself  for  some  distance  until  it  approaches  the  asymptotic 
state  where  the  self- similar  character  prevails  again.  This  state  may 
be  reached  roughly  when  the  sublayer  swallows  up  the  mass  flow  con¬ 
tained  in  the  initial  velocity  profile  immediately  after  the  corner.  Thus, 
the  complicated  transient  phenomenon  of  the  overall  viscous  layer  is 
now  replaced  by  the  changing  proportion  of  the  sublayer  and  the  shear 
layer;  the  latter  is  absorbed  in  the  former.  In  this  argument,  a  sim¬ 
plifying  assumption  has  been  made  that  the  velocity  profile  of  the  outer 
shear  layer  does  not  change  appreciably  for  the.  range  considered.  In 
other  words,  additional  mass  entrained  within  the  shear  layer  and  the 
displacement  effect  of  the  sublayer  on  the  shear  layer  are  neglected. 
Further  simplifications  may  be  made  by  approximating  linearLy  the 
average  initial  profiles  immediately  after  the  corner.  The  sublayer 
will  be  first  assumed  fully  turbulent,  from  the  corner.  Since  it  is 
assumed  to  preserve  its  self-similar  character  all  the  way,  (i/7)th 
power  law  is  adopted  for  the  velocity  profile.  Correspondingly,  the 
heat  transfer  rate  can  be  calculated  with  the  flat  plate  equations.  The 
proper  outer  edge  conditions  correspond  to  those,  of  the  initial  profiles 
obtained  right  after  expansion.  These  outer  edge  conditions  will  pro¬ 
vide  the  asymptotic  solution.  Having  now  established  the  turbulent 
asymptotic  solution  to  the  problem  at  some  distance  away  from  the 
corner',  and  also  having  established  the  laminar  region  immediately 
downstream  of  the  corner,  it  is  now  necessary  to  connect,  the  two 
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solutions.  In  this,  a  study  of  transition,  is  necessary  in  order  to  obtain 
the  complete  solution. 

The  transitional  region  is  analyzed  following  references  18  and 
19.  Use  is  made  of  the  same  simplifying  assumption  that,  the  profiles 
of  the  outer  shear  layer  do  not  change  appreciably  for  the  range  con¬ 
sidered.  The  growth  of  the.  sublayer  in  the  transformed  plane,  is  deter¬ 
mined  by  solving  the  momentum  equation  using  t.he  empirical  skin  friction 

Ve 

expression  in  terms  of  Re.  =  rr: — 7—  . 

8*  u-  =$ 


dx. 


(32) 


(33) 


Correspondingly,  the  thickness  of  the  transitional  sublayer  is  deter¬ 
mined  by  the  following  relation: 


lam 


(34) 


6* 

where  the  constants  D  and  G  are  obtained  by  matching  c,  and  (-5 — V raxas 
with  those  of  the  laminar  flow  at  the  assumed  transition  point,  Eq.  (32) 
is  integrated  numerically,  yielding  the  distribution  of  the  momentum 
thickness.  From  the  local  thickness  of  the  sublayer  and  from  the 
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the.  corner,  were  also  included  in  reference  2.  In  order  to  make  the 
boundary  Layer  turbulent  ahead  of  the  corner,  trips  were,  placed  at 
the  nose  portion  of  the  modeL  in  order  to  induce  transition.  The  pres¬ 
sure  distribution  over  the  mod.el  was  measured  with  the  trips  and 
indicated  that  there  was  no  change  from  the  one  obtained  without  the 
trips.  Therefore,  in  these  analyses  the  pressure  distribution  presented 
in  reference  2  was  used  in  calculating  the  heat  transfer  before  and  after 
the  corner.  The.  heat  transfer  measurements  are  presented  in  Fig.  4 
in  the  form  of  Nn//Re  versus  x  in  inches.  In  order  to  insure  that  the 
boundary  layer  before  the  corner  was  turbulent,  standard  heat  transfer 
analysis  was  used  to  predict  the  heat  transfer  before  the  corner.  The 
boundary  Layer  profiles  ahead  of  the  corner  have  been  calculated  with  the 
technique  described  in  Section  lH  and  are  presented  in.  Fig.  4a.  The 
Reynolds  number  based  on  the  momentum  thickness  calculated  for  the 
profiles  is  Re  a  208.  The.  profiles  obtained  after  the  in.visc.id  expansion 
around  the  corner  are  presented  in.  Figure  4b.  The  expansion  in  this 
case  has  a  two-fold  effect  on  Reg.  First  it  i  finds  to  increase  the 
momentum  thickness  after  the  expansion,  and  secondly  it  reduces  the 
density  due  to  the  expansion.  The  Reynolds  number  based  on  the 
momentum  thickness  after  the  expansion  is  143.  This  is  sufficiently  low 
to  enable  a  new  laminar  boundary  layer  to  prevail  for  quite  a  large  dis¬ 
tance  after  the  corner.  The  heat  transfer  after  the  corner  was  there¬ 
fore  calculated  with  the  aid  of  the  laminar  boundary  Layer  Eqs.  (17)  and 
(19)  using  the  turbulent  shear  profiles  presented  in  Fig.  4b.  The 
results  are  presented  in  Fig.  4,  and  it  is  clearly  seen  that  the  experimental 
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case,  Rea  was  calculated  to  be  750.  The  shear  Layer  profiles  obtained 
after  the  inviscid  expansion  around  the  corner  are  presented  in  Fig.  5b. 
For  this  condition  Re0  was  calculated  to  be  770.  Therefore,  in  this 
case,  the  expansion  around  the  corner  was  not  sufficient  to  make  up  for 
the  increase  in  the  momentum  thickness,  and  therefore  for  this  con¬ 
dition  it  would  be  expected  that  the  second  transition  to  turbulent  flow 


would  occur  in  a  rather  short  distance  after  the  corner. 


Since  a  laminar  boundary  Layer  does  prevail  after  the  discon¬ 
tinuity,  the  heat  transfer  has  been  calculated  according  to  reference  2 
and  is  presented  in  Fig.  6.  On  the  same  figure,  the  asymptotic  solu¬ 
tion  of  the  turbulent  boundary  layer  is  also  presented.  It  is  clearly 
indicated  that  transition  occurs  between  x  =  0.5->1.0  inches.  In  this  case 
transition  was  chosen  to  occur  at  an  Re  «  1100  (approximate  flat  plate 
value)  corresponding  to  x  -0.7S  .t.ches .  Wr.h  this,  the  transitional  heat 
transfer  was  calculated  with  *  he  a.d  of  Eq.  (32)  and  (.33),  and  is  included 
in  Fig.  b. 

The  previous  experiment ,  ar.c  the  ,  on  side  rations  presented  in 
this  paper  bear  o_.t  well  that  a  .ammai  -ublaver  js  followed  by  a  second 
transition  tc  turbulent  flow. 
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TABLE  1 


n 

fl 

2 

fl' 

2 

fi" 

_ i _ 

0 

0 

0 

0.7950 

0.2 

0.0159 

0.1590 

0,7945 

0.4 

0.0636 

0.3176 

0.7910 

0.6 

0.1429 

0.4750 

0.7818 

0.8 

0,2534 

0.6298 

0.7646 

1.0 

0.3945 

0.7802 

0.7381 

1.2 

0.5651 

0.9244 

0.7029 

1.4 

0.7638 

1.0610 

0.6615 

1.6 

0,9889 

1.1890 

0.6185 

1.8 

1.2388 

1.3087 

0.5799 

2.0 

1.5119 

1 .4216 

0.5521 

2.2 

1.8072 

1,5306 

0.5407 

2.4 

2.1241 

1.6392 

0.5486 

2.6 

2.4631 

1.7513 

0,5760 

2.8 

2.8252 

1.8707 

0,6197 

8.0 

3.2120 

2.0000 

0,6746 

3.2 

3.6259 

2.1418 

0.7344 

3.4 

4.0692 

2.2937 

0.7931 

3.6 

4.5441 

2.4577 

0.8462 

3.8 

5.0529 

2.6315 

0.8908 

4.0 

5.5973 

2.8134 

0.9260 

4.2 

6,1787 

3.0014 

0.9521 

4.4 

6,7982 

3.1937 

0.9703 

4.6 

7.4564 

3.3891 

0.9825 

4.8 

8.1539 

3.5864 

0.9901 

5.0 

8.8911 

3.7849 

0.9946 

5.2- 

9.6680 

3.9841 

0.9972 

5.4 

10.4847 

4,1838 

0.9987 

5.6 

11 .3415 

4,3836 

0.9994 

5.8 

12.2382 

4,5835 

0.9997 

6.0 

13.1749 

4.7834 

0.9999 
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V 
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0 

0 

0 

1  .1564 
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0,0232 

0.2312 

1.1564 

0^4 

0.0926 

0.4626 

1.1564 

0.6 

0.2082 

0.6938 

1.1568 

0.8 

0.3700 

0.9254 

1.1582 

1.0 

0,5784 

1,1574 

1.1634 

1.2 

0,8332 

1.3912 

1.1768 

1.4 

1.1350 

1.6292 

1.2056 

1.6 

1.4854 

1,8752 

1.2598 

1,8 

1.8860 

2.1354 

1.3506 

2.0 

2.3410 

2.4186 

1.4898 

2.2 

2.8558 

2.7352 

1.6866 

2.4 

3.4382 

3.0974 

1.9464 

2.6 

4.0986 

3.5180 

2.2682 

2,8 

4.8500 

4.0086 

2.6462 

3.0 

5.7074 

4.5794 

3.0692 

3.2 

6.6876 

5.2382 

3.5234 

3.4 

7.8090 

3.9900 

3.9952 

3.6 

9.0900 

6.8368 

4.4722 

3.8 

10.5500 

7.7786 

4.9448 

4.0 

12.2076 

8.8140 

5.4076 

4.2 

14.0816 

9.9408 

5.8576 

4.4 

16.1898 

11.1562 

6.2952 

4.6 

18.5498 

12.4580 

6,7214 

4.8 

21,1786 

13.3442 

7.1390 

3.0 

24.0930 

15,3132 

7.5500 

3.2 

27.3094 

16.8640 

7.9566 

5.4 

30.8440 

18.4956 

8.3602 

5,6 

34.7130 

20.2080 

8.7622 

5.8 

38.9326 

22.0004 

9.1630 

6.0 

43.5186 

23.8732 

9.5636 

37 


FIG.  I  SCHEME  FOR  ANALYZING  CORNER ,  TURBULENT  AHEAD  OF  CORNER. 
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FIG.  3  LAMINAR  HEAT  TRANSFER  DISTRIBUTION 
FOR  A  15*  SHARP  CONE- CYLINDER. 
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FIG. 4b  TURBULENT  BOUNDARY  LAYER  PROFILES 
AFTER  EXPANSION  FOR  A  24°  BLUNTED 
CONE- CYLINDER. 
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SIT  BOUNDARY  LAYER  PROFILES 


EXPANSION  FOR  A  15®  BLUNTED 
LINDER . 
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FIG.  5b  TURBULENT  BOUNDARY  LAYER  PROFILES 
AFTER  EXPANSION  FOR  A  15°  BLUNTED 
CONE -CYLINDER. 
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NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  3. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 
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Upon  comparing  Eqs,  (9)  and  (10)  there  results  a  relation  So~^0 


Furthermore, 


Vi  =  7T^~T  [hon_(ho'hw)uj^’  ki=  '  7T"  S 

2  H  U  61  O  O  W  2  U  6l 


U"  +  f  Si'  -  f  'i  l  +  2f  "i  1  =  o 

3-  O  ?  O  a  O  * 


£i  (0)  =:  £1,{0)  =  o  ,  f  i'(-n)  -  *  -  B  as  co 
2  2  2 


H*  *  W  -  V«i  =  0 


gi  (0)  =  0  ,  gi{r|)  **  T|  -  3  as  r|-»oo 
2  2 


R  =  Z° 

ueJ> 


fi"'  +  f  fV'  -  2f  'V  +  3f  "fj  =  0 
o  o  u 


fi<0)  =  fi'(O)  =  0 


fj'tn)  -  ns  -  2PH  +  0a  +  1  as  T1  - 


6 


5 

Also  included  i.e  Van.  Driest's  flat  pLate  solution  in  which 
assumptions  (1)  and  (2)  are  relaxed.  Correction  to  these  effects  can 
be  easily  made.  However,  as  is  evident  from  Fig.  2,  present  results 
show  satisfactory  agreement  with  the  experiment  for  most  engineering 
purposes. 

In  order  to  further  substantiate  the  applicability  of  this  theory, 
the  results  of  the  analy  sis  are  compared  with  the  experimental  results 
of  reference  20. 

The  pertinent  flow  quantities  for  the  15°  sharp  cone-cylinder 

are: 


Stagnation  pressure 

Stagnation  temperature 

"Wa LI -to -  stagnation  temperature  ratio 

Re 

s 

Free  stream  Mach  number 

Mach  number  outside  the  boundary 
Layer 


=  139  psia 
=  71 0°R 

-  refer  to  Fig.  3 
=  3.74x  10B 
=  5 

=  3.8  before  the  corner 
5.12  after  the:  corner 


From  these  data,  Eq.  (19)  finally  results  in 


.  ioa  =  iiHil  +  1.693  +  •••  (20a) 

Re  a  a/x" 

s 


where  pi  is  assumed  zero  and  x  is  measured  in  inches,  and  the  initiaL 
profiles  are  linearized.  The  results  are  plotted  in  Fig.  3.  Again,  the 
results  show  satisfactory  agreement  with  the  experiment.  The  experi¬ 
mental  dataievte  been  replotted  in  terms  of  the  Nusselt  number  and 
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ReynoLds  number  defined  in  this  paper. 

Finally,  consideration  is  given  to  the  applicability  of  the  sub¬ 
layer  model.  The  concept  of  the  sublayer  model  is  entirely  dependent 
on  the  characteristics  of  the  initial  velocity  profile  at  the  corner 
which,  after  being  expanded  abruptly,  has  a  slipping  portion  on  the 
wall  right  at  the  corner.  Hence,  the  corner  becomes  a  singular  point 
at  which  the  sublayer  may  be  assumed  to  start. 

If,  on  the  other  hand,  this  slipping  portion  is  missing  in  the 
initial  profile  the  sublayer  cannot  be  expected  to  start  even  if  other 
flow  quantities  may  change  discontinuously  at  that  point.  Such  an 
example  may  be  given  of  the  flow  field  immediately  downstream  of  a 
localized  injection  over  a  flat  plate.  It  can  be  shown  that  even  if  one 
starts  with  the  sublayer  model,  one  finally  arrives  at  a  merged  layer 
of  Couette  type  flow  which  indicates  that  the  sublayer  does  not  start 
at  the  end  of  the  injection  region. 


SECTION  III 
TURBULENT  FLOW 

For  turbulent  flows,  mathematical  treatment  becomes 
increasingly  difficult.  Nevertheless,  it  will  be  shown  here  that  the 
sublayer  model  coupled  with  some  simple  physical  arguments  results 
in  a  sufficiently  accurate  prediction  of  the  heat  transfer  rate.  For 
this  case  the  boundary  layer  ahead  of  the  corner  will  be  turbulent.  Two 
cases  will  be  considered  here.  The  first  case  will  be  for  the  condition 
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has  an  effect  of  reducing  the  turbulence  level  in  the  boundary  layer. 

This  is  similar  to  the  reduction  in  turbulence  level  obtained  by  the 

1  5 

contraction  section  of  a  wind  tunnel.  Morkovin  has  made  a  detailed 
study  of  the  effect  of  a  sudden  expansion  on  a  turbulent  boundary 
layer.  In  this  study,  the  expansion  was  produced  on  a  flat  plate  by 
means  of  a  wedge.  Detailed  profile  measurements  of  pitot,  pressure 
and  velocity  fluctuations  before  and  after  the  expansion  were  made.  The 
results  indicated  that  the  assumption  of  isentropy  in  the  boundary  layer 
during  the  expansion  is  valid.  The  measurements  of  the  velocity  fluc¬ 
tuations  in  percent  of  local  mean  velocity  before  and  after  the  corner 
indicated  a  decrease  close  to  the  walL,  and  an  increase  towards  the  out¬ 
side  part  of  the  boundary  layer.  Therefore,  the  results  of  references  1, 

2,  and  15  justify  the  inviscid  expansion  assumed  in  this  part  of  the 
paper.  In  each  case  the  momentum  thickness  after  the  expansion  is 
calculated  from  the  density  and  veLocity  profiLes  obtained  after  the 
expansion  with  the  following  equation: 

9  =  J  7%" (1*  Tr)dri  <31) 

o  Fe  e  e 

C.  Analysis  of  Boundary  Layer  Downstream  of  Corner 

Reference  1  has  indicated  that  the  effect  of  a  large  pressure 
expansion  around  the  corner  enables  the  growth  of  a  new  laminar  boundary 
layer  starting  at  the  beginning  of  the  cylinders.  It  is  also  believed*'  that 
the  extent  of  the  new  laminar  boundary  layer,  and  the  transition  to 
turbulent  flow,  is  controlled  by  the  high  turbu lence  level  present  in  the 
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data  is  welL  represented  by  the  theory  presented  here.  On  the  same 
figure,  the  turbulent  solution  for  the  heat  transfer  for  the  asymptotic 
state  is  also  presented  for  verifying  that  the  boundary  layer  obtained 
after  the  expansion  is  laminar.  It  is  expected  that  this  new  laminar 
boundary  layer  will  prevail  for  quite  a  distance  downstream  of  the 
corner  since  the  Reynolds  number  based  on  the  momentum  thickness 
is  low. 

Case  2,  (Reference  6) 

Stagnation  pressure  =  57.3  psia 

Stagnation  temperature  =  660°R 

Wall-to- stagnation  temperature  =  0.924 

ratio 

Re  =  5.13  x  10E 

s 

Free  stream  Mach  number  =  3.04 

Mach  number  outside  the 

the  boundary  layer  =  1.88  before  the  corner 

2.39  after  the  corner 

The  results  of  these  tests  have  been  taken  from  reference  6. 
The  body  tested  consisted  of  a  spherically  capped  15°  half  angle 
cone  with  a  bluntness  ratio  of  Ro/r=0,6  (Fig,  6). 

For  the  test  conditions  presented  above,  the  Reynolds  number 
was  sufficiently  high  so  that  a  turbulent  boundary  layer  was  well 
established  ahead  of  the  corner.  In  order  to  substantiate  this,  the 
turbulent  heat  transfer  has  been  calculated  ahead  of  the  corner,  and 
is  presented  in  Fig,  6  with  the  measurements.  It  is  clearly  seen  that 
the  experimental  results  are  well  presented  by  the  turbulent  solution. 
Fig,  5a  presents  the  boundary  layer  profiles  before  the  corner.  In  this 
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FIG.  6  HEAT  TRANSFER  DISTRIBUTION  FOR  A  15°  BLUNTED 
CONE- CYLINDER. 


